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ABSTRACT 

A lithium depletion boundary is detected in HIP 112312 (GJ 871.1 A and B), a ~ 12 Myr old pre-main sequence 
binary system. A strong (EWw 300 mA) Li 6708 A absorption feature is seen at the secondary (~M4.5) while no 
Li 6708 A feature is detected from the primary (~M4). The physical companionship of the two stars is confirmed 
from common proper motions. Current theoretical pre-main sequence evolutionary models cannot simultaneously 
match the observed colors, brightnesses, and Li depletion patterns of this binary system. At the age upper limit of 
20 Myr, contemporary theoretical evolutionary models predict too slow Li depletion. If true Li depletion is a faster 
process than predicted by theoretical models, ages of open clusters (Pleiades, a Persei, and IC 2391) estimated 
from the Li depletion boundary method are all overestimated. Because of the importance of the open cluster age 
scale, development of self-consistent theoretical models to match the HIP 112312 data is desirable. 

Subject headings: (stars:) binaries: visual — stars: pre-main-sequence—stars: abundances — stars: individual 
(WW PsA, TX PsA) — open clusters and associations: individual (/3 Pictoris moving group, 
Pleiades, a Persei, IC 2391) 


1. INTRODUCTION 

Gravitationally contracting low mass (M < 1M 0 ) stars are 
fully convective (Hayashi track evolution) until they near the 
main sequence. As the contraction proceeds, the stellar tem¬ 
perature rises and when the core temperature reaches ~ 2.5 x 
10 6 K, depletion of lithium begins via 1 Liip.af He reaction 
(Ventura et al. 1998, and references therein). Since the lithium 
burning process is very sensitive to the core temperature (burn¬ 
ing rate oc T 20 . Bildsten et al. 1997), once Li burning starts, 
it will consume all available lithium in a very short timescale 
(< 10 Myr). However, more massive stars (M > 1M 0 ) form 
radiative cores in their early evolution that results in a separa¬ 
tion of the surface layers from the hotter inner regions. As long 
as the temperature at the base of the convective layer is below 
the Li burning temperature, surface Li can be preserved. On 
the other hand, if stellar mass is too low (< 0.065M 0 ), a star 
(in fact, a brown dwarf) can never attain a high enough core 
temperature to burn Li; hence such brown dwarfs should al¬ 
ways show strong Li absorption at 6708 A. For the higher mass 
brown dwarfs (mass range of 0.065Mg -0.072Mg), lithium 
eventually can burn (Basri 2000). 

Because of the fast Li burning process, the location of the 
Li burning boundary in a coeval stellar group should be well 
defined. The Li depletion boundary (LDB) of older clusters 
stays at 0.065Mg, since that is the mass below which brown 
dwarfs never get hot enough to burn lithium. This implies that 
the LDB age determination only works up to an age of about 
250 Myr which is the age when the interior of a ~ 0.065M 0 
brown dwarf just creeps above the lithium burning temperature. 
LDBs of younger clusters will be located at higher masses than 
the LDB of old clusters because lesser mass stars take more 
time to reach the lithium burning temperature by gravitational 
contraction than more massive stars do. Once the luminosity 
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of cluster stars at the LDB is known, the LDB location can be 
compared to that of theoretical models from which we can es¬ 
timate the cluster age. To date, LDBs have been detected in 
three open clusters — Pleiades (Basri et al. 1996; Stauffer et al. 
1998), a Persei (Stauffer et al. 1999), and IC 2391 (Barrado 
y Navascues et al. 1999). Interestingly, ages determined from 
the LDB method for these clusters are all ~ 50% larger than 
traditional upper main sequence fitting ages (Table 1), and the 
lithium ages have generally been adopted as correct. If so, then 
traditional interior models for the upper main sequence must be 
revised. 

We have detected the LDB in a pre-main sequence (PMS) 
binary system HIP 112312 from medium and high resolu¬ 
tion spectroscopic observations. To date, this is the only sys¬ 
tem other than open clusters whose LDB is detected. The 
HIP 112312 binary system has a well determined distance from 
Hipparcos (d = 24 pc) and BVI C JHK S photometric magnitudes. 
Using these data, we can check if theoretical models of photo¬ 
metric isochrone ages and Li depletion ages are self-consistent. 
In this paper, we compare frequently used theoretical models 
with parameters derived from our observations. 

2. DATA 

2.1. Observed Data 

We obtained spectra of WW PsA (=GJ 871.1 A) and TX PsA 
(=GJ 871.1 B) as part of an extensive on-going search for 
young and nearby stars to Earth (Song, Bessell, & Zuckerman, 
ApJ in prep.). High resolution spectra were obtained using an 
echelle spectrograph on the Nasmyth-B focus and medium res¬ 
olution spectra were obtained with Double Beam Spectrograph 
(DBS) on the Nasmyth-A focus of the Australian National Uni¬ 
versity’s 2.3 m telescope. The red channel of the DBS cov¬ 
ered the spectral range 6500-7450 A at a measured resolu- 
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Table 1 

Li depletion boundaries in open clusters 


cluster age (Myr) Li depletion boundary Reference 


name 

UMS n 

LDB 

sp. 

M,,. 

(R-I)c 


IC 2391 

35 

53 

M5.0 

10.25 

1.91 

Barrado y Navascues et al. (1999) 

a Persei 

70 

90 

M6.0 

11.47 

2.12 

Stauffer et al. (1999) 

Pleiades^ 

85 

125 

M6.5 

12.20 

2.20 

Stauffer et al. (1998) 


“ These are upper main sequence fitting ages 

t LDB of the Pleiades was first discovered by Basri et al. (1996) at slightly earlier spectral type 


tion of 1.2 A (0.55 A/pixel). DBS spectra, displayed in Fig¬ 
ure 1, have ~ 5000 counts per pixel in the vicinity of 6700 A. 
From these spectra, one can clearly see the non-detection and 
detection of the Li 6708 A absorption feature from the pri¬ 
mary and secondary, respectively, and strong Ha emission from 
both stars. Eight orders of the echelle covered portions of the 
spectra between 5800 and 7230 A. At orders containing the 
Ha and Li 6708 A lines, the measured resolution was 0.45 A 
(0.17 A/pixel). 

All spectra were reduced following standard procedure using 
IRAF. We also obtained BVI C photometric magnitudes of the bi¬ 
nary as part of a larger photometry program (Shobbrook et al., 
in preparation) for young stars identified in our on-going spec¬ 
troscopy program. The HIP 112312 components have JHK S 
magnitudes from the 2MASS 2nd release database. Spectro¬ 
scopic and photometric data are summarized in Table 2. 

2.2. Derived Data 
2.2.1. T eff and spectral types 

Some colors of the HIP 112312 binary from BVI C JHK S are 
slightly inconsistent compared to normal M dwarfs. For exam¬ 
ple, the B — V of the two stars are bluer by about 0.1 mag than 
normal M dwarfs. The lower gravity of a young M dwarf com¬ 
pared to the mean gravities of ZAMS M dwarfs could account 
for much of this difference and a few stars in the Taurus-Auriga 
region studied by Kenyon & Hartmann (1995) have similarly 
blue B — V colors for their V—I colors. The Pleiades M dwarfs 
are also bluer in B — V than one would expect for their V -I or 
V-K colors (Stauffer et al. 2002, in prep.). This is probably 


Table 2 

Characteristics of HIP 112312 binary system 



GJ 871.1A 
WW PsA 

GJ 871.1 B 
TX PsA 

RA (2000) 

22:44:57.94 

22 : 45 : 00.04 

DEC (2000) 

-33: 15:01.7 

-33 : 15 : 26.0 

Sp. type 

M4e 

M4.5e 

V 

12.16±0.03 

13.42±0.02 

B-V 

1.50±0.02 

1.58±0.02 

V-I c 

2.78±0.02 

3.04±0.01 

J 

7.78±0.02 

8.68±0.02 

H 

7.14±0.04 

8.05±0.04 

K s 

6.91±0.04 

7.77±0.04 

R.V. (km/sec) 

+2.5±0.6 

-1.7±2.4 

EW(Ha) A 

-6.4 

-4.3 

EW(Li) mA 

< 30 

290 

p a (mas/yr) 

179.8±5.2 

179.8±2.9 

ps (mas/yr) 

-129.2±3.8 

-126.0±3.0 

Ti05 index 

0.38 

0.34 


just a feature of young K and M dwarfs. The V -1 and 
V-K colors are quite consistent and correspond to spectral 
types around M4 and M4.5 for HIP 112312 A and B, respec¬ 
tively. 

We also measured the Ti05 spectral index (= F\1 126 — 
7135]/E[7042-7046]) for HIP 112312 A and B together with 
that for GJ 644 (an M3 standard) following Reid et al. (1995). 
The measured Ti05 index for GJ 644 (TiO5=0.48) was in good 
agreement with its standard value and with its V — I color. 
Based on the Ti05 and spectral type relation of Reid et al. 
(Spectraltypesubclass = —10.775 x Ti05 + 8.2), the Ti05 in¬ 
dices for HIP 112312 A & B again indicate M4.1 and M4.5 
spectral types, respectively. From colors and Ti05 indices, 
we assign M4 and M4.5 spectral types to the primary and sec¬ 
ondary, respectively. 

The effective temperatures of HIP 112312 A and B from their 
V —I and V-K colors using an empirical temperature calibra¬ 
tion and model colors (Bessell 1991) are estimated to be 3150 K 
and 3030 K respectively, and their uncertainties are not larger 
than 100 K. 

2.2.2. Lithium Abundance 

Estimation of Li abundances from observed Li 6708 A equiv¬ 
alent width requires information on Li curves of growth. From 
synthetic model spectra (Allard et al. 2001 and private commu¬ 
nication with P. Hauschildt for updates) for 7)^=3000, 3100, 
& 3200 K and logg=4.0 & 4.5 with a range of Li abundances 
(logN(Li) = —99, 0.05, 0.1, 0.3, 0.5, 0.75, 1.0, & 1.3), we con¬ 
structed a set of Li 6708A curves of growth (Figure 2). We note 
that the indicated curves of growth take into account the effect 
of strong TiO band absorption features around 6708A which 
can alter the lithium equivalent width measurement depend¬ 
ing on the location of the pseudo-continuum level (Pavlenko 
2002, priv. comm.). When our Li curves of growth are com¬ 
pared with Pavlenko’s (2002, priv. comm.) model data which 
also takes the effect of TiO bands into account, our curves 
of growth differ by < 15 mA at common Li abundance points 
(log N(Li) = 1.0 & 1.3). Using the estimated effective tempera¬ 
ture of the secondary from section 2.2.1 (3030 K) and the mea¬ 
sured equivalent width of 290 mA, from Figure 2 we conserva¬ 
tively estimate the Li abundance of the secondary to lie in the 
range log N(Li) = 0.85- 1.10. If we adopt log n(Li) = 3.31 as 
the interstellar Li abundance, then HIP 112312 A and B cur¬ 
rently have only < 0.01 % and 0.35-0.62 % of their initial Li 
contents. 

2.2.3. Luminosities 

It is quite straightforward to calculate luminosities of the 
HIP 112312 binary system from Hipparcos measured parallax 
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FIG. 1.— Part of the DBS spectra of the HIP 112312 system. The dotted line 
represents the primary spectrum and the solid line is for the secondary. The 
secondary spectrum shows steeper slope than the primary spectrum not be¬ 
cause of the temperature difference, but presumably due to stronger molecular 
bands (i.e., CaH). 


(42.4 ±3.4 mas) and measured photometric magnitudes. In cal¬ 
culation of bolometric magnitudes, we use m—M = 1.86 for the 
distance modulus, M hoi = 4.75 for the Sun, bolometric correc¬ 
tion of 2.73 and 2.77 magnitude at K band for the primary and 
the secondary, respectively (Bessell et al. 1998). The calculated 
bolometric magnitudes are 7.78 & 8.68 and the calculated lu¬ 
minosities of A & B components are log (L/Lq) = -1.22 and 
-1.58, respectively. Uncertainty in luminosity is estimated to 
be A log (L/Lq) = 0.160 which is mainly due to the Hipparcos 
distance uncertainty. 


2.2.4. Space Motions 

Physical companionship of the HIP 112312 components is 
verified from the common proper motions of the primary and 
the secondary (Table 2). The proper motions listed in Ta¬ 
ble 2 are weighted means of SPM 2.0 (Platais et al. 1998) and 
UCAC1 (Zacharias et al. 2000) catalog values. Proper mo¬ 
tions of the primary from Hipparcos catalog (p a = 183.12 ± 
2.50 mas/yr, p,g = -118.87 ±2.21 mas/yr) are in fair agreement 
with the Table 2 values. 

From our echelle spectra, we determined weighted mean ra¬ 
dial velocities of the primary and secondary from a line fitting 
method (see Table 2). We also estimated the rotational veloc¬ 
ity (vsin i) of the primary to be (~ 40 km/sec) but were un¬ 
able to do so for the secondary because of its low S/N spec¬ 
trum. Using the measured radial velocities, Hipparcos dis¬ 
tances, and proper motions, we derived galactic space mo¬ 
tions (U, V, W ) with respect to the Sun of HIP 112312 by using 
the formulation of Johnson & Soderblom (1987); ( U,V,W) = 
(—11.8 ± 1.2,—19.0 ± 1.6,—11.0±0.9 )km/sec for the primary 
and (-13.8±1.5,-19.1 ± 1.6,-7.2 ±2.2)km/sec for the sec¬ 
ondary. Based on the similar space motions, proximity to Earth 
(24 pc), and very young ages (< 20 Myr, see Figure 3 and 
discussion in section 2.3), the HIP 112312 binary components 
may be members of the (3 Pictoris moving group (Zuckerman 
et al. 2001). The [3 Pictoris group was defined as nearby stars 
that have UVW velocity each within ~ 2 km/sec of those of 
/3 Pictoris, (-10.8,-16.4,-8.9), and age consistent with 12 Myr 
(Zuckerman et al. 2001). 



log N(Li) 

Fig. 2.—Li 6708 A curves of growth for logg = 4.3 constructed from a set of 
synthetic model spectra (Allard et al. 2001 and see text for details). The rect¬ 
angular box represents an acceptable range of Li abundance for the secondary 
star of HIP 112312. 


2.3. Discussion and Comparison to Models 

We compare pre-main sequence (PMS) stellar evolutionary 
models with our observed HIP 112312 binary data by plotting 
HIP 112312 A and B on a log T e ff versus log (L/Lq) plot along 
with three PMS models (Figure 3). The high luminosities of 
A and B clearly indicate their young age (< 20 Myr); however 
the scatter among the estimated ages from different models is 
rather large. 

The three theoretical models plotted in Figure 3 also predict 
Li depletion patterns for PMS stars. Since we know how much 
Li currently remains in HIP 112312 A and B components (sec¬ 
tion 2.2.2), we can directly compare the observed Li contents 
to the predicted values (Table 3). 

Zuckerman et al. (2001) identified ~ 30 co-moving stars near 
Earth, the (3 Pictoris moving group. By plotting late-type mem¬ 
bers on a color magnitude diagram with theoretical isochrones, 
they estimated an age of the group as ~ 12 Myr. 



Fig. 3.— Hertzsprung-Russell diagram of the A and B components of 
HIP 112312 with isochrones from three different PMS theoretical stellar evo¬ 
lutionary models (Baraffe etal. 1998; D’Antona & Mazzitelli 1997; Siess etal. 
2000). From top to bottom, 10, 20, & 100 Myr isochrones are plotted for each 
model. The dotted luminosity error bars indicate the case if either A or B 
should itself happen to be an equal mass binary. However, to the Adaptive 
Optics resolution at Keck, neither the primary nor the secondary is a binary 
(Macintosh et al., in prep.). If HIP 112312 A & B are ~ 12 Myr old, then 
their masses are ~ 0.32 Mq and ~ 0.18M©, respectively, based on the above 
models. 
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Table 3 

Predicted and observed Li contents for stars with logL/L 0 =-1.22 (primary) and logL/L 0 =-1.58 (secondary). 


Model 

age 

1 Li content (% 

of initial content) 


(Myr) 

Primary 

Secondary 

Observed values 

<20 

<0.01 

0.35-0.62 

Baraffe et al. (1998) 

10 

~ 98 

100 


20 

~ 0.7 

~ 56 

D’Antona & Mazzitelli (1997) 

10 

~ 85 

100 


20 

<0.01 

~ 5 

Siess et al. (2000) 

10 

~ 98 

100 


20 

<0.01 

~ 69 


Note. — Among different theoretical models, D’Antona & Mazzitelli (1997) fits the observed data best and all LDB ages of open clusters were based on the 
D’Antona & Mazzitelli model. However, even at the age upper limit of the HIP 112312 binary system, the D’Antona & Mazzitelli model still predicts too slow Li 
depletion (~ 5 % of Li content for the secondary corresponding to > 600 mA of Li 6708 A equivalent width). 


Recently, Ortega et al. (2002) estimated a kinematic age of 
the (3 Pictoris moving group by tracing positions of individual 
members back in time. They found that the (3 Pictoris mem¬ 
bers occupied the smallest volume about 11.5 Myr ago; hence 
they assigned a kinematic age of 11.5 Myr. We obtain some¬ 
what similar results when we trace positions of the currently 
known members back in time using constant velocity trajecto¬ 
ries (i.e., no galactic potential treatment. Song et al. in prep.). 
Many members were more concentrated into a central clus¬ 
ter about 12 Myr ago although the overall size of the moving 
group has not been reduced as much as claimed in Ortega et al. 
(2002). The reasonably good agreement between photometric 
isochrone age and dynamical age of the (3 Pictoris group in¬ 
dicates that ages estimated from theoretical isochrones can be 
trusted in general for ~ 10 Myr old stars. 

Based on the facts that the HIP 112312 A component lo¬ 
cates on the Hertzsprung-Russell diagram above the 10 Myr 
isochrone and that HIP 112312 AB may belong to a large 
~ 12 Myr old (3 Pictoris moving group (Zuckerman et al. 2001; 
Ortega et al. 2002), we believe that the true age of HIP 112312 
cannot be older than 20 Myr. At the age upper limit of 20 Myr, 
all three theoretical models predict too slow Li depletion (here 
we assume that for given ages, theoretical model effective tem¬ 
perature and luminosity calculations are more reliable than 
lithium depletion calculations). If the true age of HIP 112312 
is younger than ~ 12 Myr (as may be hinted from Figure 3), 
then inconsistency between photometric data and Li depletion 
timescales of contemporary theoretical models becomes more 
severe. If the true age of HIP 112312 is not significantly dif¬ 
ferent from ~ 12 Myr, thus requiring a faster Li depletion pro¬ 
cess than those predicted by contemporary models to fit the ob¬ 
served data, then ages of open clusters (Pleiades, a Persei, and 
IC 2391) estimated from the LDB method are all overestimated. 
Since an age scale of open clusters is critical and tied into many 
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areas of astronomy, developing self-consistent theoretical mod¬ 
els to match the HIP 112312 binary system data is an urgent 
task. 

Lastly, it is worthwhile mentioning that a low mass sec¬ 
ondary component (M4.5) of a previously identified /3 Pic¬ 
toris member (V343 Nor, K0V) also shows a similarly strong 
Li 6708 A feature (Song et al, in prep.). V343 Nor B and 
HIP 112312 B have almost the same spectral types based on 
the Ti05 spectral indices and photometric colors. Study of 
V343 Nor B corroborates the idea of HIP 112312 being a mem¬ 
ber of the (3 Pictoris moving group and the location of the LDB 
must be at ~M4.5 for 10-20 Myr old stars. 

3. SUMMARY AND DISCUSSION 

We have identified the Li depletion boundary in a < 20 Myr 
old pre-main sequence binary system, HIP 112312 A & B. 
From its galactic space motion, young age, and proximity to 
Earth, the HIP 112312 binary system may be a member of the 
~ 12 Myr old (3 Pictoris moving group. At the conservative 
age upper limit of 20 Myr, all contemporary theoretical evolu¬ 
tionary models predict too slow Li depletion. Incorrect model 
Li depletion rates may cause overestimated open cluster ages 
(Pleiades, a Persei, and IC 2391) from the lithium depletion 
boundary method. Because of the importance of the open clus¬ 
ter age scale, development of self-consistent theoretical models 
to satisfy HIP 112312 binary system data is desirable. 
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